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ABSTRACT 

We estimated photospheric velocities by separately applying the 
Fourier Local Correlation Tracking (FLCT) and Differential Affine Veloc- 
ity Estimator (DAVE) methods to 2708 co-registered pairs of SOHO/MDI 
magnetograms, with nominal 96-minute cadence and ~ 2" pixels, from 46 
active regions (ARs) from 1996-1998 over the time interval t^^ when each 
AR was within 45° of disk center. For each magnetogram pair, we com- 
puted the reprojected, average estimated radial magnetic field, 5/?; and 
each tracking method produced an independently estimated flow field, 
u. We then quantitatively characterized these magnetic and flow fields 
by computing several extensive and intensive properties of each; exten- 
sive properties scale with AR size, while intensive properties do not de- 
pend directly on AR size. Intensive flow properties included moments of 
speeds, horizontal divergences, and radial curls; extensive flow properties 
included sums of these properties over each AR, and a crude proxy for 
the ideal Poynting flux, Sr = \u\B]^. Several quantities derived from 
Bji were also computed, including: $, total unsigned flux; R, a measure 
of the unsigned flux near strong- field polarity inversion lines (SPILs); and 
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Yl^R- Next, using correlation and discriminant analysis, we investigated 
the associations between these properties and flares from the GOES flare 
catalog, when averaged over both and shorter time windows, of 6 and 
24 hours. Our AR sample included both flaring and flare-quiet ARs; the 
latter did not flare above GOES Cl.O level during Among magnetic 
properties, we found R to be most strongly associated with flare flux. 
Among extensive flow properties, the proxy Poynting flux, Sr, was most 
strongly associated with flare flux, at a level comparable to that of R. All 
intensive flow properties studied were more poorly associated with flare 
flux than these extensive properties. Past flare activity was also associ- 
ated with future flare occurrence. The largest coefficients of determination 
from correlations with flare flux that we performed are ~ 0.25, implying 
no single variable that we considered can explain variations in average 
flare flux. 

Subject headings: Sun: atmospheric motions; Sun: magnetic fields; Sun: 
flares 



1. Introduction 



(A version of this paper with much higher-resolution images is available at 



http://solarmuri.ssl.berkeley.edu/~welsch/public/manuscripts/welschjetaL20090904.pdf) 



Solar flares are intermittent releases of energy from the solar corona, in the form 
of accelerated electrons and ions, and strongly enhanced electromagnetic radiation 
over a wide range of wavelengths, from radio to X-ray. Consideration of the grav- 
itational, thermal, and magnetic energy densities in the corona suggests that only 
the coronal magnetic field contai ns sufficient energy to account for the estimated en- 
ergy release on flare time scales (jForbesll2000l ). Observations of stable, non-potential 
(e.g., twisted and/or sheared) magnetic field structures associated with coronal en- 
er gy releases (e.g., sigm oids associated with coronal mass ejections [CMEs] studied 
by lCanfield et al.lll999l ). have led to wide ac ceptance of the storage and release model 



of solar flares (see, e.g., [Linker et al.ll200ll ). In this model, magnetic free energy — 
energy above that of the potential magnetic field matching the same normal-field 
boundary condition at the base of the corona — is built up and stored in the corona 
on relatively long time scales (from several hours to weeks) and then suddenly re- 
leased on much shorter time scales (minutes to a few hours). This naturally poses 
two questions: (1) What processes lead to the buildup of magnetic free energy in the 
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corona? (2) And what processes lead to the sudden release of free magnetic energy 
in flares? 

Free magnetic energy in the corona can increase by: (a) flows that carry new flux 
from beneath the photosph ere into the c orona, and/or (b) photospheric flows that 
twist, shear, and/or cancel (IWelschI |2006| ) already-emerged fields. Obser vations and 
theore tical models have demonstrated the viability of these mechanisms. iLeka et al. 
(119961 ) convincingly demonstrated that twisted (non-potential) fields emerge through 
the photosp here, carrying free magnet i c energy w ith them. Simulations of this emer- 
gence (e.g.. iMagara fc Longcopell200ll : lFanll200ll ) produced sig moidal (S- or iii verse- 



S-shaped) field lines that may become visible in X-rays (e.g.. iFan et al.ll2003l ). and 



resemble the C ME- associate d sigmoids studied by ICanfield et al.l (119991 ). In already- 
emerged fiel ds. iBrown et al.l (120031) found rotating sunspots in some eruptive active 



regions, and IWelsch et al.l (120041 ) observed shear fiows along the eastern neutral line 



in NOAA AR 8210 prior to a fiare and CME. In addition, rotating, shearing, and 
cancelling p hotospheric fiows have been used to iniect free magnetic energy in simu- 



lations fe.g..lLvnch et a 



2002; Tokman fc Bellan 



20081 
2OO2I : 



Roussev et al 



2004 



Antiochos et al 



regio n fiows have been conducted (e.g.. 



Keil et al 



Amari et al. 2003; Lionello et al 



19991). While case studies of active 



19941 : iLi et al.l I2004J : lOeng et al 



20061 ). investigations of the quantitative properties of photospheric fiows and their 
relationships to fiares in many active regions are lacking. 

To remedy this shortcoming, we employed tracking techniques to estimate fiows 
from magnetogram sequences of a few dozen active regions (ARs) that exhibited 
varying levels of fiaring activity. We then investigated correlations between properties 
of the estimated fiows and the occurrence of solar fiares. We also studied relationships 
between fiare occurrence and several properties of AR magnetic fields. Of particular 
note, two magnetic field properties that have been previously related to solar fiares 
can serve as benchmarks for fiow-fiar e associations: the tot al unsigned magnetic fiux 
$ of the AR at the photosphere (e.g.. lLeka fc Earned 120071 ): and Schrijver's (2007) R 



value, a measure of the unsigned fiux near strong-field polarity inversion lines (SPILs). 

As an aside, we note that Schrijver's R value is very similar to a measure de- 
veloped earlier by Falconer et al. (2003, 2006) — the length of SPILs, which should 
be strongly correlated with R by virtue of the strong- field threshold. In fact, Barnes 
et al. (2009) report a correlation coefficient of 0.95 between these two quantities for 
their AR sample. Falconer et al. (2003, 2006) were among the first to quantitatively 
characterize the strong correlation between the presence of strongly non-potential 
photospheric magnetic fields and strong gradients in both radial and line-of-sight 
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(LOS) mag netic fields . Whil e the work Falconer et al. (2003, 2006) dealt with CMEs, 
and that of ISchriiverl (120071) d ealt with flares, these two phenomena are clearly re- 



lated (see, e.g.. [Andrews 



20031 ) . As a practical matter, we found Schrijver's measure 



simpler to compute than that of Falconer et al. (2003, 2006). 

We are aware of studies of helicity evolution in active regions which suggest that 
the flux o f magnetic helicity ac ross the photosphere might play a role in solar flares 
(see, e.g., LaBonte et all 2007 ). Tests of flow estimation methods using simulated 



data (IWelsch et al.l 120071 ). however, demonstrated that helicity estimates are often 
error-prone. Accordingly, we plan to investigate flare-helicity relationships in a future 
study, in which possible sources of error can be carefully examined. 

In the next section, we describe our data set of magnetogram sequences from 
a few dozen ARs, and in §3], we describe how we estimated velocities from these 
sequences. In §H we describe how we quantitatively characterized the intensive and 
extensive properties of the magnetic fields and estimated flows. In ^ we explore 
associations of these magnetic and flow properties with flares. Finally, in ^ we 
discuss what we have learned from our study. 



2. Data 

We chose to estimate flows from magnetic evolution observed in sequences of 



MDI (jScherrer et al.lll995l ) full-disk (~ 2" pixels), 96-minute nominal cadence, LOS 
magnetograms. Compared to available magnetogram sequences with higher resolution 
and higher cadence (including some in which the full magnetic vector is estimated), 
the near-ubiquitous coverage of these full-disk LOS magnetograms is a tremendous 
advantage for tracking. After we began this study, the SOI/MDI team recalibrated 
field strengths in the MDI full-disk, 96-minute database; as our analysis had already 
progressed substantially, we chose to stick with the Level 1.8.0 data, which only 
included the plate scale correction^ 

While sequences of white-light data could have been tracked instead, it is unclear 
precisely how white-light evolution is related to the underlying flow field. Further, 
evolution of the photospheric magnetic field is presumably more directly related to 
evolution of the coronal magnetic field — which powers flares — than evolution of 
photospheric white-light intensity. In addition, the relatively short lifetime of white- 



http://soi.stanford.edu/magnetic/Levl. 8/ 
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light intensity features — most notably granules — places stringent limits on data 
cadence. Since magnetic structures appear much more persistent than white-light 
features, tracking with slower cadences is possible (though not necessarily desirable, 
a point discussed further below). 

We selected several dozen ARs and AR remnants from MDI synoptic maps from 
1996 to 1998, from solar minimum into the ascending phase of the most recent cycle. 
Our sample was unbiased with respect to flare and CME productivity, but not with 
respect to AR morphology: we chose ARs that were, subjectively, both roughly bipo- 
lar and isolated (at the photosphere). Some of the AR remnants that we selected and 
tracked lacked sunspots, and therefore also lacked a NOAA AR designation; these 
have been excluded from the results presented here. The AR sample we use here is 
comprised of Nar = 46 ARs. 

Properties of our AR data set are presented in Table [H We only tracked each 
AR while it was within about 45° of disk center as it transited the disk; the start 
and stop times for tracked intervals are shown in the table. To present a quantitative 
measure of the variability of flare activity within our sample, we have also included a 
measure of each AR's average daily GOES flux from flares at or above the Cl.O level, 

^45 = (1005^^) + 105(*') + 1.05(^))/^45 , (1) 

where T45 is the time interval (in days) that the AR was tracked across the portion 
of the disk within 45° of disk center, an d >5'^^-' is the s um of GOES flare significands 



in the i-th GOES class over r45, after lAbramenkd (120051 ). The units of are 
fi W m~^ day~^. To associate flares with ARs in our sample, we used flare source 
designa tions from the NO AA flare catalogj^ although this catalog is known to contain 



biases (jWheatlandll200ll ). Flares associated with our ARs, during the time interval 
that each AR was tracked, include: zero A-class flares, 119 B-class flares, 154 C- 
class flares, 15 M-class flares, and two X-class flares. While the true frequency of 
occurrence of A- and B-class flares should exceed that of C-c lass flares, based upon 



the distribution of flare peak fluxes (e.g., ICrosby et al.l 119931 ). weak flares are often 



not ascribed to any AR in the GOES catalog, and we have only included flares with 
a source AR in our study. Of our 46 ARs, only 24 flared above the GOES Cl.O level 
while within 45° of disk center; of the remaining 22, a further 8 were identified as the 
source region of at least one sub- C-class flare. 



^http: / /www.ngdc.noaa.gov/stp/SOLAR/ftpsolarflares.html=i^xray 
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3. Methods 

The MDI instrument measures the LOS field strength, -Blos? averaged over each 
pixel. To estimate the radial field, Br, from the observed -Blos, cosine corrections 
were applied to the LOS field in each pixel, 

5r = 5los/cos(7) , (2) 

where 7 is the heliocentric angle from disk center to each pixel, equal to the inverse 
sine of the distance (in pixels) from disk center to each pixel divided by the solar 
radius (in pixels). 

To compensate for foreshortening, triangulation was used to interpolate the Bji 
data — regularly gridded in the plane of the sky, but irregularly gridded in spher- 
ical coordinates (6*, 0) on the solar surface — onto points (x, y) corresponding to a 
regularly gridded Mercator projection of the spherical surface. This projection was 
adopted because it is conformal (locally shape-preserving), which is necessary to en- 
sure displacements measured by our tracking algorithm were not biased in direction. 
The price of conformality is distortion of scale, such that apparent displacements are 
exaggerated by a factor of the secant of apparent latitude; but this is easily corrected 
after tracking. Details of our deprojection procedure are discussed in Appendix \M 



We separately used the Fourier Local Correlation Tracking (FLCT: lFisher fc Welsch 



20081) and the Differential Affine Velocity Estimator (DAVE; ISchuckl 120061 ) codes to 
produce two estimates of flows for each pair of successive magnetograms in our se- 
quences. Prior to tracking, each pair of magnetograms was cross-correlated and then 
co-registered to sub-pixel accuracy using cubic-convolution interpolation. For FLCT, 
pixels below a 20 G threshold in |-B_r|, where B^i is the average magnetic fleld for each 
magnetogram pair, were not tracked, nor were pixels farther than 7 = 45° from disk 
center. DAVE estimated flows for many pixels that fell outside of these thresholds, but 
a fllter was applied to the raw DAVE results that zeroed out flows from pixels where 
FLCT did not estimate flows. Because the physical length corresponding to Mercator 
pixels shrinks toward the poles, all estimated flows were corrected for the distortion 
in magnitude introduced by the Mercator projection. Moving cropping windows were 
applied to the deprojected magnetograms to extract subimages containing the ARs 
as they crossed the disk. 

While the magnetogram cadence was nominally At = 96 minutes, data gaps 
did exist. Flows estimated from magnetogram pairs with At > six hours apart were 
excluded from the analyses of flow properties described in ^ With At = 96 minutes , 



pixels of As ~ 1400 km, and flows of ~ 0.3 km s ^ (cf., iHathaway et al.l |2002| ) 
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displacements would be ~ 1.2 pixels. We expect that this cadence effectively averages 
away ffows with characteristic time- and length-scales shorter than a couple of hours 
and a couple of megameters. This implies that the flows we measure ought to be 
long-lived and large-scale relative to gr anular flows typically trac ked from sequences 
of intensity images in white-light (e.g.. [November fc Simonlll988l ) or G-band. 



Demoulin fc Bergen (120031 ) asserted that tracking methods applied to magne- 



togram sequences estimate an apparent flow u, with 



(3) 



and where v is the plasma velocity, r is the normal to the magnetogram surface, 
and h denotes components in the magnetogram surface. Here, we are using Bpt as 
an estimate of Br, the true normal (radial) field. Substitution of equation into 
the normal component of the magnetic induction equation produces a continuity-like 
equation for the normal componen t of of the ma gr ietic field and the flux tran sport ve- 
locity u JPemoulin fc Bergei][2003l : ISchuckll2005h . | Demoulin fc Bergeii (120031 ) pointeci 
out that knowledge of u and B is sufficient to compute the fluxes of magnetic en- 
ergy and helicity across the magnetogr am surface. Howeve r, recent tests of tracking 
methods on MHD simulation data by IWelsch et al.l (120071 ) cast doubt on the accu- 
racy of methods that exclu sively track th e normal component of the magnetic field. 
A foUowup investigation by ISchuckl (120081 ) using the same MHD simulation data sug- 
gest that some tracking methods estimate v/j much more accurately than u — they 
produce a biased estimate of the horizontal plasma velocity. This supports the view 
that line-of-sigh t tracking methods are i nsensitive to the vertical flows that drive 
flux emergence (IDemoulin fc PariatI l2009l ) . However, these methods may produce a 
reasonable estimate of the fluxes as sociated with f ootpo int t wisting in matu re active 
regions. Tests conducted by both IWelsch et al.l (120071 ) and ISchuckl (120081 ) suggest 
that velocity estimates that employ the full magnetic vector — not just the LOS data 
that we use here — would be more accurate. Currently, however, sequences of vector 
magnetograms are too rare for statistical studies relating magnetic evolution to flares. 

For tracking each active region's data cube, both FLCT and the DAVE require 
selecting a spatial windowing parameter, a. In principle, this choice can be made 
objectively, by selecting the value of a that maximizes agreement between the ob- 
served magn etic evolution ABn/At b etween magnetogram pairs, and the expected 
relationship (IDemoulin fc Bergen l2003l ) between this evolution and the observed flow 
u, 

Vh ■ (ufiif) = -ABR/At , (4) 
where Bji is the average of Bji from each magnetogram in the pair. (But LCT schemes 
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are not formally consistent with this continuity- like equation; see Schuck 2005.) In 
practice, the cadence in our study was so low that equation (jl]) was relatively poorly 
obeyed (compared to tests with other data sets with higher cadence) for any choice 
of a. Tests that quantified agreement of DAVE results with equation (jl]) using slope- 
fitting, as well as linear (Pearson) and rank-order correlations, suggested that a in 
the range of 8 — 10 pixels maximized consistency across these measures. Accordingly, 
cr = 9 pixels was used with the DAVE, and cr = 8 pixels was used with FLCT. We 
expect this windowing leads to further spatial smoothing of flows, averag ing over flows 



with c haracteristic length scales less than several megameters. Following IWelsch et al 



(I2OO4J ). positive- and negative- flux pixels were tracked separately with FLCT, with 



these separate estimates combined into all-polarity flow maps. 



The FLCT code used here was slightly modified from that described by lFisher fc Welsch 



(2008). The modified versioiifl removes the mean intensity from the subimages prior 



to cross correlating them. A Gaussian spatial frequency filter was used in the FLCT 
runs, with a roll-off wavenumber in x of /Cx.roU-off = 0.25 m.?i:dk,x)i a nd similarly for y. 



This has been found to improve accuracy in tests with known shifts (iFisher fc Welsch 



20081 ). and to ameliorate FLCT's tendency to underestimate speeds. 



Our full data set consists of A^aow = 2708 sets of flows estimated by both FLCT 
and DAVE. A typical FLCT flow field is shown in Figure [TJ This example was taken 
from AR 8038, an acti ye region which produced three flares and CMEs as it crossed 



the disk in May 1997 (iLi et al.ll2009l ). The time shown on the plot corresponds to 
the end of the 96-minute interval over which the flows were estimated. DAVE flows 
from tracking the same data are shown in Figure [2l While there are broad regions 
of agreement between the directions of FLCT and DAVE flows, there are also some 
regions where the flows are in the opposite direction. Our focus here is not to compare 
FLCT and DAVE flows, but rather to test the sensitivity of flow-flare correlations (or 
lack thereof) to the method chosen to estimate flows. Nonetheless, in Figure [3l we 
show scatter plots relating the x- and |/-components of FLCT and DAVE flows from 
the same time interval, to illustrate how the flows compare in this case. We only 
plotted results for the set of pixels in which both methods estimated flows. On the 
same plot, we also show the linear (Pearson) and rank-order correlations between the 
flow estimates. Clearly, the methods' flow estimates are significantly correlated, but 
do not agree in detail. For the ~ 16.5 million pixels for which both FLCT and DAVE 
estimated velocities, the linear and rank-order correlations of the x-components of the 



^available at 



http://solarmuri.ssl.berkeley.edu/^fisher/public/software/FLCT/C-VERSIONS/ 
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flows were 0.67, while for the y-components, these correlations were 0.60 and 0.61, 
respectively. Also, the linear and rank-order correlations between FLCT and DAVE 
speeds, when both are weighted by the unsigned field strengths, are higher still — 
0.85 and 0.76, respectively. This suggests the flows tend to agree more closely in 
strong field regions. 

Comparing speeds shows a systematic difference. Over the same set of ~ 16.5 
million pixels, FLCT's mean and median flow speeds were 0.069 and 0.062 km sec~^, 
while for DAVE, mean and median flow speeds were 0.112 and 0.099 km sec~^. In 



tests with synthetic magnetograms from ANMHD simulations, IWelsch et al.l (120071 ) 
found that an older version of FLCT than that used here underestimated speeds by 
about ~ 30%, which would suggest a DAVE/FLCT slope of 0.7"^ ~ 1.4. Here, the 
DAVE/FLCT speed ratio is nearer 1.6. 

These flow speeds are significantly less than flows estimated by other methods, 
e.g., supergranu l ar flow speeds determined from center-to-limb Doppler variation by 



Hathawav et al.l (120021) . They are, however, only slightly slower than flows estimated 



by lChae et al.l (120041 ) for 96- minute cadence data; our slower speeds might result from 



our larger spatial windowing, ~ 17" (see their Fig. 7). 



4. Quantifying Magnetic Field and Flovi^ Properties 

4.1. Typical Fields and Changes 

To familiarize the reader with our data set, we present some statistics regarding 
typical magnetic fields and field changes in the approximately 21 million pixels in 
which \Bji\ exceeded 20 G (of about 85 million total pixels in our data set). Over this 
set of pixels, the median and average cosine-corrected, unsigned, pixel-averaged field 
strengths {Bui were 63 G and 113 G, respectively. The median and average unsigned 
rates of change, | Ai?/j|/At were 3 x 10^^ G s^^ and 4 x 10~^ G s~^, respectively, which 
correspond to unsigned changes of 17 G and 23 G, respectively, for the nominal 
At = 5760 s time interval between magnetograms. Perhaps unsurprisingly, this is 
roughly the quoted noise level for MDI magnetograms]^ The linear and rank-order 
correlation coefficients between unsigned \Bii\ and \ABii/At\ were 0.36 and 0.41, 
respectively — hence, absolute changes tend to be larger in stronger fields. The linear 
and rank-order correlation coefficients for \Bji\ and \{ABf(/At)/Bf(\ were -0.30 and 



http://soi.stanford.edu/magnetic/Levl. 8/ 
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-0.39, respecively, implying relative changes tend to be smaller in stronger fields. The 
rank-order correlation between signed Bji and AB^i/At was 0.05, which amounts to 
a small but statistically significant bias: negative fields were slightly more likely to 
get more negative, and positive fields were slightly more likely to get more positive. 
This imbalance, however, cannot be sustained indefinitely. Our ARs all have NOAA 
designations (and so possess sunspots), and therefore contain relatively newly emerged 
flux; perhaps older flux in decaying active regions behaves differently. 



4.2. Field Properties 

In this and the foHowing subsections, we describe aspects of magnetic and flow 
flelds that we calculate to relate to flare activity. We distinguish between "extensive" 
properties, which depend directly on AR size (e.g., integrated quantities), and "inten- 
sive" properties, which do not depend directly on AR size (e.g., averaged quantities). 
This is in analogy to thermodynamics, in which extensive quantities scale with the 
system size (e.g., total thermal energy, or particle number), while intensive proper- 
ties do not (e.g., temperature, or particle density). As a rule, quantities that would 
increase if an AR's length scales were doubled (all else held equal) are extensive, e.g., 
total unsigned flux; quantities that would not increase are intensive, e.g., average 
magnetic fleld strength. Although our characterization of $ as extensive is based 



upon a hypothetical doubling of AR area, we note that iFisher et al.l (119981 ) found 



that $ is in fac t strongly positively correlated with (and nearly a linear function of) 



AR area. While iFisher et al.l (119981 ) used the term "extrinsic" for $ and other integral 



fleld properties, we prefer extensive for its more direct analogy to thermodynamics. 

While our analysis centers on properties of estimated flows, we also characterized 
properties of the time-averaged estimated radial magnetic fleld Bn corresponding to 
each pair of fleld estimates used for tracking, as well as properties of changes in the 
magnetic fleld that do not require tracking. These magnetic fleld and fleld evolution 
properties are described below, and then summarized in Table |2j 

In the subsequent statistical analyses of magnetic fleld properties, it will of- 
ten be appropriate to use magnetic fleld values scaled in a way to compensate for 
the latitude-dependent area distortion introduced by the Mercator projection. This 
rescaling is necessary because the physical area represented by Mercator-projected 
pixels shrinks toward the poles (the effect which causes Greenland and Antarctica to 
appear unrealistically large in Mercator projections of the Earth's surface); details 
are given in Appendix |X1 Without rescaling, fleld strengths in Mercator pixels that 
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represent a smaller area on the Sun would be given equal weighting in our analyses. 
The apparent area da of a Mercator-projected pixel — equal to the MDI disk-center 
pixel area, (~ 1.4 Mm)^ — is related to Mercator-corrected pixel area da via 

dd{6) = da cos^ 9 , (5) 

where 9 denotes the latitude of each pixel (cf., co-latitude 9 often used in spherical 
coordinates). Pixel values used in calculations involving local magnetic properties or 
local image structure (e.g., the field strength threshold for tracking, or contours of 
Br) do not require Mercator correction. In contrast, pixel values in global quantities 
(e.g., sums, averages, and higher moments of quantities) must be weighted by cos^ 6*. 
Because we frequently refer to Mercator-corrected magnetic field values, we define 

I3r = Br cos^ 9 , (6) 

and 

Bl^Blcos'9. (7) 



Leka fc BarnesI (120071 ) demonstrated that the total unsigned flux was strongly 
associated with flaring, so we calculated this extensive quantity from each average 
estimated radial field, 

<^ = J2\BR\dd = Y,\BR\da, (8) 

where the sum runs over the cropping window. 

For each magnetogram, we also computed the first four moments of the estimated 
radial magn etic field, B^: mean (i ?/?), variance a'^[BR], skewness C[Br], and kurtosis 



k[Br], after Leka fc BarnesI ( 2003 ). We denote this family of moments Ai[BR], 



M[Bf 



(Br) 


mean 


ct'IBr] 


variance 


([Br] 


skewness 


4Br] 


kurtosis 



(9) 



Since each of these moments does not, of necessity, scale with active region size, these 
descriptors of AR magnetic fields are intensive. 

Because fiare processes might be nonlinearly related to photospheric field strength 
(because, for instance, magnetic energy density is proportional to magnetic field 
squared), we also computed 
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where, as above, the sum runs over all pixels in each average magnetogram. While 
B]^ has units of magnetic energy density, we do not ascribe any particular physical 
meaning to the sum of this quantity; it is included in our analysis to provide baseline 
comparisons for other quantities. (We also investigated sums of higher powers of B^, 
but linear and rank-order correlations of summed B]^ with summed l-B/j]^ and Bj^ are 
all above 0.9, implying these variables contained little additional information.) 



Schrijven (120071 ) correlated the total unsigned flux, R, near strong-field PILs 



(SPILs) with flaring, so we also computed R for each magnetogram in our data set. 
This entails: (1) creating two separate positive and negative bitmaps, representing 
regions of strong average field {\Br\ > 150 G) in each polarity; (2) computing the 
product of dilated versions of these bitmaps, to find loci with strong, opposite po- 
larity fields in close proximity; (3) computing a SPIL- weighting function, Wspil, by 
convolving the product bitmap with a normalized Gaussian of 15 Mm FWHM; (4) 
summing the unsigned flux in the SPIL-weighted magnetogram. 



R 



spil\Br 



(11) 



A few of these steps are illustrated explicitly in IWelsch fc Lil (120081 ) . Some magne- 
tograms in our study have i? = 0. While R quantifies one aspect of the internal 
arrangement of flux within an AR — how much lies near strong-field polarity inver- 
sion lines — it also scales with active region size, albeit with the length L of each 
AR, not area A. Accordingly, we also characterize R as extensive. 

In addition to $ and i?, we calculated the time derivatives $ and R between 
averaged magnetograms (each in [Mx sec~^]), to investigate the correlation between 
changes in these quantities and flares. If the flux emergence rate per unit area, in 
Mx sec~^ cm~^ , were uniform over the Sun, then th e nearly linear dep endence of $ 



with the area A over which magnetic flux is present (IFisher et al.lll998l ) would imply 
a nearly linear relationship between $ and $ — essentially, the larger an AR's $ 
i s, the larger the cross-section for emergence into that AR. In fact, iHarvey fc Zwaan 
( 119931 ) showed that emergence rates are enhanced over the background level in areas 
where flux has already emerged. Hence, we characterize $ as extensive, and expect 
$ and $ to be correlated. To account for this, we also calculated intensive measures 
of the rates of change of $ and R, the flux-normalized $/$ and R/ R, to compensate 
for correlations between $ and and R and R. 



Wand (120061 ) reported changes in the center-of-flux positions of each polarity 



associated with five fiares (X- or M- class) in (5-spot regions. Accordingly, we also 
computed the center-of-flux (COF) positions (Mercator corrected) of each AR polarity 
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in each average cropped magnetogram, 



R\ 



(12) 



± 



where the ± sums run over pixels of the corresponding polarity. We then computed 
the mean polarity separation vector, 



dx 



X 



(13) 



which points from x+ to x_, and found the rate of change (dx) (in [km sec~^]) of 
dx between pairs of average magnetograms. We then computed the "center-of-flux 
divergence," 

Dc:oF = (dx)-dx/|dx| , (14) 



and normalized the change by the initial dx, as well as the "center-of-fiux shear," 



^coF = |i--((dx)xdx)|/|dx| 



(15) 



implied by these changes. The unit vector r is normal to the projected magnetogram 
plane. Both Dcof and Scof have units of velocity. A negative (resp., positive) Dqof 
corresponds to centers of flux growing closer (resp., f arther) along the liii e separating 



them. Based upon flux emergence simulations (e.g., lAbbett et al.ll2000l ). one might 
expect Dcof > with the emergence of new flux, and, conversely, Dcof < with 
flux cancell ation. If cente rs of flux have a tendency to move closer together as a result 
of flares, as IWang (120061 ) found, then flares should be anti-correlated with Dcof- A 
non-zero Scof corresponds to displacements of centers of flux in opposite directions 
perpendicular to dx. 



4.3. Flow Properties 

Here, we describe how we quantifled apects of our estimated flows; these prop- 
erties are also summarized in Table [2l For each pair of magnetograms we track, we 
typically have two independent flow estimates, by FLCT and DAVE, for each of sev- 
eral thousand pixels. From these estimates, we computed several intensive measures 
of flow properties. As with Bji, we computed the moments of each flow fleld, A^[m]. 
Because photospheric flows in regions of stronger 1 5^ | are expected to impact coronal 
evolution more strongly than regions with weak |-B_r|, we also computed normalized 
and unnormalized flux- weighted moments, A^[u|(i$|/$„] and A^[M|i?ij|], respectively. 
We note that u\Bji has units of an electric fleld, is related to the horizontal electric 
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field. In the normalized expression, \d^\ represents the (Mercator-corrected) flux in 
each pixel, and is the total unsigned flux only from those pixels for which each 
method (FLCT or DAVE) estimated flows. (Averaging over the entire 2D arrays for 
each velocity estimate — which contain many zero velocities [e.g., from weak-fleld 
regions] could bias the moments if, for instance, low-flux ARs have fractionally more 
weak-fleld pixels.) We also computed several related extensive quantities: summed 
speeds and fleld- weighted speeds, and ^wl-BijI, respectively; and the sum of a 
kinetic-energy like quantity, X]""^) fleld-weighted version, J2'^'^\^r\, respec- 

tively. A grayscale map of u'^\Bii\ is shown in the upper-left panel of Figure HI 

We also computed the sum of a quantity that, dimensionally, has units of an 
energy flux, 

(16) 

with u in cm sec~^ (our default in this paper is km sec~^) and Bj^ in G^, the summand 
has units ergs cm~^ sec~^. Multiplying by MDI's disk-center pixel area da, ~ 2 x 10^^ 
cm^, then yields units of power, ergs sec~^. This sum represents an extensive measure 
of the transport of magnetic energy at the photosphere. The true radial Poynting 
flux Sr across the photosphere and into the corona can be derived from an estimated 
flow fleld u and a vector magnetogram (jDemoulin fc Berger 2003 ) via. 



Sr = (l/47r)^[(B,,-B,>,-(v;,-B;,)5,] (17) 

= -(l/47r)5^(u-B,)B, , (18) 

which depends sensitively on the relative orientation of the flow and magnetic flelds. 
Our simple proxy Sr completely ignores this depend ence. As long as m agnetic flelds 



are not exactly vertical, and observations reported by lFisher et al.l (119981 ) suggest they 
tend not to be, then strong horizontal flelds will tend to occur where there are strong 
vertical flelds. Further, flux emergence and submergence occurs on PILs; but our 
proxy is small along PILs, where the radial fleld vanishes. But flux from emerging 
(submerging) features should diverge (converge), which can be detected (see §4.4^ 
below). Consequently, our proxy, crude though it may be, could be statistically 
related to actual energy fluxes. A grayscale map of Sr is shown in the upper-right 
panel of Figure HI 

To characterize flow properties arising from the vector nature of our estimated 
flows, we computed several additional intensive properties of the flow. For each flow 
fleld u, we computed the absolute and signed horizontal divergence, and the unsigned 
radial component of the vorticity (curl), and then computed moments of these dis- 
tributions. When computing derivatives for these quantities, Mercator corrections 
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were applied to account for varying pixel scales. As above, horizontal refers to direc- 
tions perpendicular to the normal (radial) direction, f . In addition, we computed the 
flux-weighted signed and unsigned horizontal divergences, unsigned radial curls, and 
moments of these distributions. We excluded derivatives in pixels with any of the 
four-nearest neighbor pixels below the absolute field-strength tracking threshold, to 
avoid contributions to divergences and curls from flow discontinuities between tracked 
and untracked pixels. Weighted divergence and curl maps are shown in the left and 
right (resp.) panels in the bottom row of Figure |H In addition, we computed ex- 
tensive versions of these quantities, totaling signed and unsigned, unweighted and 
field-weighted divergences and curls. 

Numerical simulati ons of coronal evolution leading to coronal mass ejections often 



use s hearing flows (e.g., lAntiochos et al.l 119991 ) or converging fl ows with cancellatio n 



le.g. 



Linker et al.ll2003l ). or some combination of the two (e.g., iRoussev et al.ll2004j ). 
along PILs to both inject energy into the coronal magnetic field and form eruptive 
structures. In addition, flux eme rgence is expected to lead to diverging flows near 



PILs (see, e.g.. lAbbett et al.ll200Cll ). Basically, these flow patterns can be distinguished 
by their relationships with the local gradient in normal fleld: shearing velocity vectors 
point along normal-fleld contours, and converging or diverging flows near PILs are 
parallel or anti-parallel to normal-fleld gradients there. Accordingly, to characterize 
the presence of shearing or converging/ diverging flows near PILs in an automated 
way, we also decomposed flow estimates {u} into components along x^, a unit vector 
that points along the local horizontal gradient of Bji, 

X, = WhBr/IWhBrI (19) 

and along Xc, a unit vector that points along contours of Br, 

^, = rxVhBR/\VhBR\ . (20) 

We denote these (signed) gradient and contour flows as {ug} and {uc}, respectively. 
Examples of gradient and contour flows are shown in the top panel of Figure 0, as 
blue and red vectors, respectively. We computed A1[mc], A1[Mg], 7M[mcM'J* and 
I (i$ !/$„]. We also summed these quantities for each flow fleld to create extensive 
measures of flows related to shearing and convergence. 

We also computed PIL-weighted contour and gradient flows, A4[WpilUcBr], 
and M.\WpilUcBr\, respectively. The PIL- weighting function used here, Wpn, was 
meant to identify only regions of magnetograms near PILs. It was computed in 
the same way as the SPIL weighting function Wspil (which was used to determine 
R; see above), but with different parameters. The polarity bitmaps used for Wpn 
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had a threshold of only 40 G, while the threshold for the bitmaps used to create 
WspiL maps was 150 G; and the overlap-map of the dilated polarity-mask bitmaps 
was convolved with a normalized Gaussian of ~ 7 Mm FWHM for WpjL, versus 
15 Mm WspiL- Weighting by the signed magnetic flux, introduces a cross-PIL 
parity, emphasizing regions with oppositely directed flows across the PIL, and de- 
emphasizing regions where flows on opposite sides of the PIL are similar. Grayscale 
maps of (WpjLd^yweighted gradient and contour flows are shown in the left and right 
(resp.) bottom panels of Figured In each of these bottom panels, contours of Wpn 
for Bp from the top panel are overplotted in red. For some magnetograms, Wpil = 
everywhere, and for some others, [Wpnu) = everywhere, i.e., regions with \Bp\ 
above the tracking threshold and non-zero PIL-map weighting do not necessarily 
overlap. We note that the product of a horizontal velocity estimate with the radial 
magnetic field corresponds to a horizontal electric field perpendicular to both; hence, 
the left and right grayscale maps effectively show the components of the horizontal 
electric field along contours and gradients (resp.) of Bp. 



4.4. Expected Flow Properties 

As an aside, we sought to confirm some expected correlations withi n the dataset 



Of no te, we expected diverging flows when flux was emerging (see, e.g., lAbbett et al 



2000l ). and converging flows when flux was cancelling. Assuming that flux emergence 
(cancellation) leads to a positive (negative) $, we computed the linear and rank-order 
correlations between FLCT's disk-passage averaged signed, flux-weighted divergence 
and disk-passage averaged $/$, and found correlation coefficients of 0.61 and 0.67, 
respectively. Using DAVE's disk-passage averaged signed, flux-weighted divergence 
in the correlations gave values of 0.73 for both linear and rank-order correlations. 
We also computed the linear and rank-order correlations between the center-of-flux 
divergence, Dcofi and and found correlation coefficients of 0.45 and 0.48, re- 

spectively. 

We also looked for evi dence that shearin g along PILs is associated with flux 



emergence, as predicted by iManchesten (120071 ). First, we correlated emerging flux 



with flux-weighted and unweighted measures of curl, but coefficients were in the 
rage 0.2-0.4 — smaller than the correlations of emergence with divergence. We also 
computed rank-order and linear correlations between Scof and $/$ of 0.31 and 
0.29, respectively; these are weaker than the association of Dcof with $/$. In 
addition, the rank-order and linear correlations between disk-passage averaged, PIL- 
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and flux- weighted contour flows (see the discussion below equation [20]) and disk- 
passage averaged $/$ were 0.13 and 0.36, respectively, for FLCT, and 0.24 and 0.42 
for DAVE, respectively. Our measures of shearing motions were generally more weakly 
correlated with flux emergence than our measures of diverging motions were. 



Tests with synthetic magnetog rams conducted bv IWelsch fc Fisherl (120081 ) with 



FLCT verified the supposition by iDemoulin fc Bergerl (120031 ) that estimated flow 



magnitudes \uc\ along contours of Bj^ would be underestimated in comparison to flow 
magnitudes \ug\ along gradients of Bfi. We did not find, however, evidence of a strong 
bias in flow speeds along and perpendicular to contours in our data set: averaging 
over the ~ 16.5 million pixels for which both FLCT and DAVE estimated flows, the 
mean and median contour speeds \uc\ were 2% - 5% higher than gradient speeds \ug\ 
for both FLCT and DAVE. Differences between real and synthetic magnetograms 
might explain the discrepancy between these results: real data are more likely to 
possess fine-scale structure and asymmetries than synthetic data, and these structural 
details make it possible for tracking methods to determine the overall flow, even along 
contours. It might still be the case, however, that our estimates are in fact biased 
with respect to contours/ gradients in Br, and that true contour speeds \uc\ exceed 
true gradient speeds \ug\ by an even greater amount. 

To understand the interplay between field strengths and estimated speeds, we 
plotted the distributions of speeds u with pixel-averaged field strengths Bji for both 
FLCT and DAVE in Figure O The Figure includes speeds from all pixels for which 
both FLCT and DAVE estimates were available. Because these several million speed 
estimates were too clustered for a regular scatter plot, we binned speeds and field 
strengths, and used log-scaled shading instead of points in high-density bins, and 
overplotted individual points in low density bins. The log-scale shading for the FLCT 
plot ranges over 2.4 - 6.2, and for the DAVE plot it ranges over 2.4 - 6.0. Clearly, 
higher speeds tend to occur in weaker fields, though a range of speeds are possible 
at any field strength. Statistically, however, correlations between speeds and field 
strengths are weak, with rank-order coeffi cents for FLCT and DAVE of 0.07 and 
-0.02, respectively. 
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5. Associating Field & Flow Properties with Flares 

5.1. Disk-Passage Flare Associations 

As a first step toward understanding the relationship between flows and flares 
(if any), we investigated the correlations between each AR's average GOES flux from 
flares of class C or greater, gross magnetogram properties (e.g., $,-R, and their time 
derivatives), and flow properties, averaged over the tracked component of the AR's 
disk passage. This coarse averaging in both space and time eliminates many detailed 
aspects of magnetic fleld structure and evolution that might be relevant to flare 
occurrence. Nonetheless, any correlations that remain despite this gross averaging 
ought to be robust. The values for flare flux used in these correlations are those listed 
in the right-most column of Table [H 

In Figures [7] and [8], we compare rank-order and linear (respectively) correlation 
coefficients between whole-disk-passage averaged magnetic fleld and flow properties, 
from both FLCT and DAVE, with whole-disk-passage averaged flare power from equa- 
tion ([T]), when averaged over the 46 active regions in our study. Quantities are labeled 
in order of decreasing distance from (0,0), the point on these plots which corresponds 
to complete lack of correlation with average flare power; hence, smaller numerical la- 
bels imply stronger correlation. Only zeroth, flrst, and second moments of magnetic 
fleld and flow properties are included in this ranking, and only the most strongly 
correlated properties are shown. (In plots that included higher moments, no mag- 
netic fleld properties were highly-ranked enough to appear. Some higher moments of 
flow properties, however, did appear; but further analyses, described below, suggest 
that these correlations were not robust.) Odd rank numbers appear below their cor- 
responding plot symbols, while even rank numbers appear above their symbols. To 
the extent that flare correlations with a given flow property from FLCT and DAVE 
flow agree more strongly, that flow property's plotted symbol is nearer to the dotted 
diagonal line in each plot (which is shown to guide the eye; it is not a flt). Quantities 
that do not require a flow estimate (e.g., R) lie on the diagonal line. We believe this 
graphical presentation of our results displays information more clearly than listing 
our results in tables would. Differences between linear and non-parametric, rank- 
order correlation coefficients could arise from lack of linearity between the correlated 
quantities, or non- normal variable distributions (e.g., numerous outliers), or our rel- 
atively low sample number Nar = 46, or some combination of these. It is notable, 
for instance, that the most highly linearly correlated quantity, R, is not among the 
most non-parametrically correlated quantities; but its "excess kurtosis" (that relative 
to the normal distribution) is 15, implying it deviates very strongly from the normal 
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distribution. 



Some of the quantities most strongly associated with our measure of disk-passage 
averaged flare flux, by either method of computing correlations, do not require esti- 
mating flows: $,-R, Yl^Ry quantities involving time derivatives of the first two. 
Some of these magnetic field - flare correlations accord w i th previous work rel ating 
flare occurrence rates with tota l flux (ILeka fc Barnes! 120071 : iBarnes fc Lekall2008l ) and 
flux near PILs (jSchriiveii 120071 : and also for CMEs, as reported by Falconer et al. 
2003, 2006) . Table B] s hows that disk-p assage averages of these quantities are corre- 
lated with each other. ISchrijverl (120071 ) has argued that large values of R are more 
predictive of the occurrence of large flares than $ alone, but our correlation analysis 
does not unequivocally support this assertion. We note that the signed, intensive rate 
of change of flux, is more strongly correlated with flare power than the unsigned 
intensive rate of change, |$|/$, implying that flux emergence is more strongly corre- 
lated with flare energy release than flux cancellation (which makes $ non-zero, but 
negative). The extensive rate of change of flux $ is more strongly correlated with 
average flare flux than the intensive rate of change of flux, which suggests that 

more flare energy is released when more photospheric flux emerges than when less 
flux does, even if the new flux is larger relative to the pre-emergence flux. 

Several extensive flow properties are significantly correlated with our measure of 
disk-passage averaged flare flux. Most notably, our proxy Poynting flux, Sji, is among 
the quantities most strongly correlated with average flare flux. 

Some intensive flow properties show relatively strong negative correlation with 
flare flux: the averages of field-weighted unsigned vorticity, \d^{r ■ V x u)|/$„, un- 
signed divergence, |(i$(V ■ u)|/<l>„. The magnetic flux - flare flux correlation is strong 
enough, however, that normalization by introduces a spurious, inverse correlation 
between these field-weighted quantities and average flare flux. We found that the un- 
normalized, field-weighted versions of these average flow properties (not shown) are 
positively correlated with flare flux (but not sufficiently so to appear on these plots 
of the most strongly correlated properties). 

Apart from these flux-normalized flow properties, most other intensive flow prop- 
erties, averaged over the few days when ARs are near disk center, are too weakly cor- 
related with average GOES flux to appear in these plots. We expect, however, that 
AR magnetic flelds are more coherent on the relatively long timescale of a few days 
compared to flows, w hich have lifetimes on the order of a few hours to a day (e.g., 
Leighton et al.lll962l ). Figures [9l and [TO] compare, for FLCT and DAVE, respectively, 
linear correlation coefficients of Br, Ux, and Uy in time — both frame-to-frame (thick 
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black, blue, and red lines, resp.) and with respect to the initial frame (thin black, 
blue, and red lines, resp.) — for ARs 7981, 8038, and 8210. Clearly, frame-to-frame 
variations are much larger for estimated flows than for B^, meaning flow proper- 
ties tend to average incoherently — toward zero — in the disk-passage averaging we 
have undertaken here. This analysis suggests flows have a characteristic decorrelation 
timescale — taken as the time required for the autocorrelation to fall to 1/e — of 
~ 6 hr. 

In Figure [9], we also overplotted flare start times, shown by letters represent- 
ing GOES flare class. (Note that the vertical placements of letters do not scale 
with actual flare fluxes.) There is no obvious correspondence between flare start 
times and field or flow autocorrelations. The differences between these three ARs' 
current-frame-to-initial-frame magnetic fleld correlations are intriguing, but investi- 
gating these variations is beyond the scope of this study. 



5.2. Flare Associations on Shorter Timescales 



5.2.1. Correlation Analysis 



The flow timescale implied by the foregoing correlation analysis — e-folding 
times of ~ 6 hr — has implications for analyzing flow-flare relationships. Accepting 
the arguments of iMcClymont fc Fisherl (119891 ). we presume that flows do not directly 
power flares, but still allow that some types of flows might tend to trigger flares. It 
is plausible that "trigger" flows have the same characteristic timescale, ~ 6 hr, as 
the general flows in our sample. Alternatively, if some types of photospheric flows 
correspond to the injection of energy (perhaps via flux emergence) into the corona 
over several hours preceding flares (thereby being the ultimate source of flare energy, 
as opposed to the proximate energy source if flows were to drive flares), then there 
ought to be a fl o w-flar e rel ationship on the corona l energy storage timescale. Work by 
Schrijver et al.l (120051 ) and iLongcope et al.l (120051 ) suggests the timescale for coronal 
relaxation via flaring and reconnection is on the order of 24 hr. Accordingly, we now 
address the relationship between flow properties and the occurrence of flares in these 
two shorter time windows, 6 and 24 hours. 



As a flrst step, we correlated magnetogram and flow properties with the average 
GOES flux from C-class and larger flares, computed via equation ([T]), but with the 
disk-passage time interval t^^ replaced with the flare window length, tq for 6 hours, 
or T24 for 24 hours. For each possible window duration, two intervals were analyzed: 
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(1) the "current" window, centered on the average time of the two magnetograms 
used to estimate flows; and (2) the "next" window of the same duration following the 
current window. For ease of reference, we label these four windows 6C, 24C, 6N and 
24N, with C and N denoting the current and next windows. To be explicit, our 6C 
window includes flares in the 6-hour interval ±3 hours from each velocity estimate, 
and 6N includes flares in the 6-hour interval 3-9 hours after that velocity estimate. 
The 24C window includes flares within ±12 hours of each velocity estimate, and 24N 
includes flares over the 24-hour interval between 12-36 after that velocity estimate. 
In space weather parlance, forecasts corresponding to the 6N and 24N windows would 
have latencies of 3 and 12 hours, respectively. Note that our time windows do not 
correspond t o the zero-latency forecast windows used by some other authors, e.g.. 



Barnes et all (120071 ) 



The scatter plots in Figures [TT] and [T2] show rank-order (R.O.) and linear corre- 
lation coefficients, for both C and N time windows, between average flare power and 
magnetogram or flow properties from FLCT (horizontal axis) and DAVE (vertical 
axis). The format of plots in these Figures is similar to that of Figures [7| and [HI Cor- 
related quantities are labeled in order of decreasing distance from (0,0), the point on 
the plots which corresponds to complete lack of correlation with average flare power. 
Odd rank numbers appear below their corresponding plot symbols, while even rank 
numbers appear above their symbols. In these plots, all moments of magnetogram 
and flow properties were ranked, and only the ten most highly ranked properties were 
shown. To the extent that flare correlations with flow properties from FLCT and 
DAVE agree, they lie along the diagonal dotted hne (not a fit). Quantities that do 
not require a flow estimate (e.g., R) lie on the diagonal line. No higher moments of 
magnetic field or flow properties beyond the zeroth, first, and second were correlated 
strongly enough with flare flux in these windows to appear in these plots. 

These plots show relationships similar to the disk-passage averaged results. As 
with the disk-passage averaged correlations, many of the quantities that are strongly 
positively correlated with average flare power do not require estimating flows — e.g., 
$, i?, ^ 5^, the average unsigned field (|-B_r|), as well as the standard deviations of 
the signed and unsigned fields, cr[-Bij] and (T[|i?i{|], respectively. Again, our proxy for 
the Poynting flux Sr is among the quantities most strongly correlated with average 
flare flux. 

The coefficient of determination is defined as r^, where r is the linear correlation 
coefficient. It quantifies how much of the variation in the dependent variable (avg. 
flare flux, here) is explained by the independent variable: when = 1, all of the 
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observed variation can be explained in terms of the correlation, and a smaller value 
of implies less explanatory power. The largest values of for the correlations with 
flare flux that we performed are ~ 0.25, implying most variation in flare flux cannot 
be explained by any single variable that we considered. 

The clustering of correlation coefficients from magnetic field and flow properties 
in Figures[TT]-[T2]hampers clear identification of which magnetic field and flow proper- 
ties are most relevant to flare activity. Hence, relationships between these properties 
and flare activity should be analyzed further. 



5.2.2. Discriminant Analysis 
We also used the disc riminant analysis (DA) code developed by Ba rnes, Leka 



and collaborators (see, e.g.. lBarnes et al.ll2007l : iLeka &: Barnesll2007l . l2003l ) to analyze 



our dataset. Briefly, given samples classified into distinct groups (in our case, time 
windows that have been classified into binary populations, either flare-associated or 
flare-quiet) and a set of observed quantities corresponding to each sample (in our 
case, magnetic field and flow properties), DA can be used to generate a discriminant 
function — a linear combination of the observed quantities — that, in a particular 
sense, best reproduces the classification of each sample into the distinct groups. This 
discriminant function can then be used as a way to predict the class to which a 
future sample will belong. While its approach of determining the linear combination 
of observed quantities that can explain variation in the samples is similar to that of 
principal component analysis, DA focuses on determining which observed quantities 
most strongly distinguish between the classes. 

In principle, with the observed quantities standardized by subtracting their 
means and dividing by their variances (using, e.g., IDL's standardize. pro), ratios 
of coefficients in the discriminant function quantify the relative discriminatory power 
of the observed quantities. Thus, DA might be used to identify the strongest dis- 
criminators, even if several of the observed quantities are cor related with each other. 



Value s of these coefficients can, however, be unstable (e.g., iHuberty fc Wisenbaker 



19921 ). implying such ratios cannot be used to rank the importance of observed quan- 
tities with coefficients of similar magnitude. Nonetheless, DA should still be able to 
distinguish between variables whose coefficients differ significantly. A variable with a 
much smaller discriminant coefficient either inherently lacks power to accurately clas- 
sify a sample, or contributes little additional discriminatory power beyond that from 
another variable in the discriminant function with which it is highly correlated. In 
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accordance with these considerations, we focus less on ranking the utihty of variables, 
versus identifying the few with strong discriminatory power in each time window. 

Among other things, DA assu mes the distribution of each observed quantity is 



normal; but iLeka fc Barnes! (120071 ) found the results of DA applied to vector magne- 
togram properties as predictors of solar flare occurrence to be relatively insensitive 
to the distributions of input properties. In addition, DA assumes that the observed 
quantities are independently measured. In our study, this is not strictly true for mag- 
netogram properties like $, R, and which were measured from the estimated 

radial field Br averaged from successive magnetograms; these quantities calculated 
from successive magnetogram pairs are not formally independent. The high frame- 
to-frame correlations for the AR magnetic fields shown in Figure [HI however, imply 
that even if formally independent magnetogram measurements had been used, the 
resulting magnetic field properties would still be highly correlated. 

We used DA to investigate the relationship between estimated magnetic field / 
flow properties and the occurrence of C-class (or stronger) flares, in each of the four 
time windows (6C, 6N, 24C, 24N) associated with each estimation of magnetic field 
and flow properties. With the nominal 96-minute cadence of estimates of magne- 
togram and flow properties in our study, several estimates might be associated with 
a single flare. Nonetheless, DA should still be able to quantify the relative power of 
observed magnetic field and fiow properties to discriminate between fiare-quiet and 
fiare-associated windows. 

While DA can identify the subset of variables with the most discriminatory power, 
the predictive power of that combined set of variables (with relative weights given by 
the discri minant function) m ust be evaluated separately. Climatological skill scores 



see, e.g., IWheatlandl 120051 ). which quantify the improvement in forecast accuracy 



over the constant f orecast derived fro m the average event frequency, can be used for 



this purpose. After IWheatlandl (120051 ) . we denote averaging over all samples by (...), 
the forecast probability by /, and the actual occurrence by x (with a value of 1 or 
for each sample, for flaring or not, respectively), with a skill score SS given by 

ss = i-{{f-xr)/{{{x)-xr) . (21) 

Perfect forecasting results in SS = 1, and scores below zero reflect worse performance 
than a uniform "climatological" prediction of flare likelihood equal to the observed 
frequency (x). While several other metrics of forecast models' predictive power exist, 
validation of a forecast model is not the object of this paper, so we opted not to 
consider additional metrics. 
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Which properties associated with each pair of magnetograms had the greatest 
discriminatory power for each time window? We note that, in addition to the mag- 
netic field and flow properties considered in the correlation analysis above, we also 
considered the average flare flux during the current time windows, 6C and 24C, as a 
property that could be included in discriminant functions for windows 6N and 24N, 
respectively. Because Sr was strongly correlated with average flare flux for many 
time windows in the analysis above, we started by using DA to construct separate, 
two-variable discriminant functions with Sr and each of the other possible magnetic 
fleld, flow, and flare (for windows 6N and 24N) variables in turn, for each of the 
four time windows. Separate versions of these discriminant functions were calculated 
using FLCT and DAVE flow properties. In this way, each property was separately 
tested "head to head" against Sr. For some windows and/ or flow estimates, it was 
found that R (and only R) had a larger discriminant coefficient than Sr. Accordingly, 
for these windows, we repeated this head-to-head analysis, constructing two-variable 
discriminant functions with R and each of the other variables in turn; and no other 
variable had a larger coefficient in the two-variable discriminant functions for these 
windows than R. These computations also identified the quantity that, for each win- 
dow, had the most additional discriminatory power: the quantity with a discriminant 
coefficient closest in magnitude to the highest. For all windows but 24C, Sr and R 
had the largest coefficients. For 24C, (|-B_r|) — the average absolute field strength, 
an instrinsic magnetic quantity — had the largest coefficient relative to Sr. 

We then constructed three-variable discriminant functions, separately for each 
additional magnetic field and flow property (also separately for FLCT and DAVE 
flows) beyond the two properties already identifled for each time window. At this 
stage, the hierarchy of discriminant function coefficients differed between the two flow 
estimation methods for all windows, though skill scores only differed by at most 0.02. 
As additional variables beyond these were considered, skill scores did not substantially 
increase, and the hierarchies of discriminant function coefficients constructed using 
variables from the two different flow estimates methods continued to disagree. Both of 
these facts suggested that including additional variables would result in diminishing 
marginal improvement in identifying magnetic fleld and flow properties associated 
with flares. Accordingly, we did not pursue this analysis further. 

Table H] summarizes the outcome of our discriminant analysis. Time windows 
are listed in the left-most column, and method used to estimate flows are hsted in 
the next column. The variables used to construct each discriminant function are 
shown in the next column. For multiple-variable discriminant functions, the discrim- 
inant coefficients of each standardized variable are listed in square brackets. These 
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results suggest that Sr and R consistently have strong discriminatory power between 
flaring and non-flaring populations, regardless of tracking method used to estimate 
5*;^. Further, given their commensurate discriminant coeflicients, our discriminant 
analysis implies that these two parameters each possess some unique predictive ca- 
pability. The (intensive) average unsigned fleld strength, (|-Bi?|), also had relatively 
large discriminatory power in many cases. 

The next four columns of Table H] show predicted and actual flare frequencies 
for each window considered. (There are many more flares in the statistics listed in 
these columns — from around a hundred to a few hundred, depending on the window 
— than are listed in ^ This discrepancy results from "overlapping measurements" : 
many average magnetic flelds and flow estimates having a common flare in their win- 
dows.) In classifying events for each window, we used the same data that were used 
to construct the discrimant function for that window, which tends to bias the success 



rate upward (lHillslll966l ). (Unbiased calculations require constructing separate dis- 
criminant functions with each sample removed, and using the resulting discriminant 
function to classify the removed sample. This is both time consuming, and, in tests 
we performed, did not signiflcantly alter event classiflcations.) The flnal column lists 
the climatological skill scores, SS, derived for the corres ponding discr i minan t func- 
tion. For comparison with the skill scores reported here, iBarnes et al.l (120071 ) used a 
large set of properties of AR magnetic flelds derived from vector magnetograms, and 
achieved an 5*5* value of 0.35 for predicting C-class or larger flares in the 24 hours 
following each magnetogram. Our differing deflnition of the future 24 hour window 
(12-hour latency in our case, no latency in theirs) precludes direct comparison of our 
results with theirs, but their results nonetheless provide useful context. 



We also note that the vector magnetograms used by lBarnes et al.l (120071 ) are prob- 
ably superior to LOS magnetogr ams for computing radial magnetic fleld strengths. 
For instance, iKnoU et al.l (120081 ) computed the skill scores for single-variable dis- 
criminant functions using total unsigned flux, computed from cosine-corrected LOS 
magnetograms at various distances from disk center; they found that the skill score 
using magnetograms within 45° of disk center (as with our data set) was a factor 
of three lower than the skill score using only magnetograms from within 10° of disk 
center. This demonstrates that using off-center LOS data to estimate Br, as we have 
done, can negatively impact skill scores. 

Skill scores were negative for the 6N window using our discriminant functions, 
implying no skill beyond that based upon event frequencies in the data set. As flares 
are infrequent, forecasting for small time windows tends to be inaccurate, as a short 
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time interval is more likely than not to be flare quiet (ILeka fc Barnesll2007l ): this might 
partly explain why forecast skill for window 24N exceeded that of 6N. For window 
24N, flare flux averaged over the current 24 hour window, ^24,current5 exhibited some 
discriminatory power, and evidently shares some predictive power with i?, as the 
latter's coefficient shrank when the former was included in the discriminant function. 
The tendency for regions that have already flare d to flare again soon afte r is referred 
to as persistence in the flare forecast literature (jZirin &: Marquettdll99ll ). 

In Table [5l we list statistical properties of the two variables for each time win- 
dow with the highest discriminant function coefficients from Table |H We first list 
the means of the raw variables for the flaring and non-flaring populations from each 
window, separately and combined; and then list the corresponding standard devia- 
tions of each variable's combined flaring and non-flaring populations. We then list 
the means of the standardized versions of each variable for the flaring and non-flaring 
populations. Results from FLCT were used for the only quantity that involves flows. 
Sr. 



6. Discussion & Conclusions 



Magnetic energy is the predominant form of energy in the low coro na, and 
free magnetic energy is thought to power flares and CMEs (IForbesI 120001 ). This 
free energy is thought to e merge into the corona in active region magnetic flelds 



(IMcClymont fc Fisherlll989). As the ph otospheric magnetic fleld should be strongly 



coupled to the plasma (e.g. . |Parkerlll984l ) . the transport of magnetic energy across the 
photosphere should be governed by equation (fT7|) : either radial (ver tical) flows trans- 
port current-carrying magnetic flux across the photosphere (e.g., iLeka et al.l Il996l : 



Abbett fc Fishei 



Brown et al 



2003 



20031). or photospheric flows act on already emerged fields (e.g.. 



Antiochos et al.lll999l ). Accurate estimates of photospheric flows. 



therefore, could be used to study physical processes that lead to the buildup of coro- 
nal magnetic energy that is released in flares and CMEs. If processes that trigger 
flares/ CMEs are related to photospheric evolution, then studying photospheric flows 
can improve our understanding of these processes. It is possible, however, that trig- 
ger processes are essentially coronal, and only related to photospheric evolution in a 
statistical sense. In the latter case, knowledge about trigger processes that could be 
gleaned from studies of photospheric evolution would be limited. 



In this study, we attempted to relate properties of magnetic flelds and flow flelds 
estimated by tracking photospheric magnetograms to the release of coronal magnetic 
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energy in the form of flares. We distinguished between extensive properties, which 
scale with AR size, and intensive properties, which do not depend directly upon AR 
size. Using correlation analysis, we identified several properties that were roughly 
equally associated with flare flux averaged over several time windows. To distinguish 
which among these properties was most strongly related to the occurrence of flares, we 
employed linear discriminant analysis. Using both methods, we found one extensive 
flow property, Sr = '^B\, to be associated with flare flux and flare occurrence at 
a level comparable to the most strongly flare-associated magnetic field property, i?, 
which quantifies the flux near strong- field p olarity iri v ersion lines (SPILs). Previous 
work by Falconer et al. (2003, 2006) and ISchrijverl (120071 ) had already linked the 
amount of flux near SPILs to CMEs and flares. 

The association of Sr with flare activity was robust regardless of whether FLCT 
or DAVE were used to estimate flows, and for a variety of time windows corresponding 
to each flow estimate. To quantify the similarity between the flows when weighted by 
-B^, we computed the linear and rank-order correlation coefficients between mflct-B^ 
and mdave-B|. for the ~ 16.5 million pixels where both estimated flows; they were 0.89 
and 0.91, respectively. We note that Sr was more strongly associated with flares than 
^ implying the square of the magnetic field strength is not solely responsible for 
association. We also found u and to be uncorrelated with each other, so the 
association with flaring does not arise from a hidden correlation between u and \Br\. 
Hence, we conclude that flow information is an essential factor in the association 
between Sr and flare activity that we have identified. We could not unambiguously 
associate any intensive property of flows with flare activity, a potentially significant 
fact we discuss in more detail below. 

The extensive or intensive scaling of flare-associated magnetic field and flow prop- 
erties is relevant to flare and CME processes. Triggers of flares and CMEs might be 
small-scale features like current sheets; for instance, globally significa r it rec onnection 
begins at current sheet s in the breakout m odel of lAntiochos et al.l (119991 ) and the 
tether cutting model of iMoore et al.l (120011 ). In contrast, the non-potential magnetic 
structures associated with flares and CMEs, such as filaments and sigmoids, are large- 
scale features. Free magnetic energy is probably stored on large scales, because the 
difference between the actual and potential magnetic fields is non-local. For example, 
the presence of a current sheet — a compact non-potential structure — can alter 
the magneti c field far from the current sheet itself. Flares occur over a range of en 
ergie s (e.g 



Hudson 



19911) . and larger ARs can produce larger flares (ISammis et al. 



2000l : iKucera et al.l 119971 ) . suggesting that peak flare energy is an extensive property 
of ARs. The effectively exponential weighting of average flare flux by flare class in 
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equation ([T]) means that the average flare fluxes we used here strongly emphasize 
the largest flares. Hence, our measure of flare activity essentially scales extensively. 
Consequently, it is perhaps unsurprising that the AR properties we found to be most 
strongly associated with flare activity are extensive. 



Figure [13] illustrates this extensive scaling. The squares plot the relationship 
between the log of the average flux from flares at or above GOES Cl.O level over the 
time interval when each AR was within 45° of disk center, versus the log of the 
average total unsigned magnetic flux in these ARs over that time. Error bars show 
RMS variance in flux for each AR over errors expected from 20 G per pixel noise in 
an individual MDI magnetogram are about the size of the plotted squares or smaller. 
Triangles show the flux in the 22 ARs that did not flare at or above GOES Cl.O level. 
(Because many A- and B-class flares are not reliably ascribed to a source AR in the 
GOES catalog, emission from these flares has been excluded.) It is clear the squares 
show a trend, which we crudely quantified by computing linear and rank-order (R.O.) 
correlations between magnetic flux and average flux from flares above the GOES Cl.O 
level, both for raw values using all ARs, and for logarithms of these values using only 
ARs that flared above this level. We also computed a least- absolute- deviation fit 
of the log of average flare flux as a function of the log of unsigned magnetic flux, 
which yielded a slope of 1.40. While our small sample size hampers our ability 
rigorously characterize this relations hip, we note that the value of the slope differs 
from the value of ~ 1.15 found by iPevtsov et al.l (120031 ) to relate total unsigned 
flux to thermal soft X-ray (SXR) luminosity. A discrepancy between the scalings 
of thermal and non-thermal (flare) SXR luminosities with AR magnetic flux might 
indicate that different physical processes generate each type of emission. Also, the 
magnetic flux - flare flux correlation might account for some of the variability in AR 
flare productivity reported in other studies re latin g flares to magnetic fi eld properties 
and/ or evolution — e.g., lAbramenkd (120051 ) and lLaBonte et al.l 120071 . An intensive 
measure of AR fiare productivity might be formed by normalizing average AR fiare 
fiux by AR magnetic fiux $ rais ed to some power a ; for a larger sample, one could 
follow the procedure adopted by iFisher et al.l (119981 ) to determine an exponent that 
removed any correlation between average flare flux and magnetic flux. Linearity in the 
magnetic flux - flare flux relationship could be consistent with uniform free energy 
input per unit flux (uniform "speciflc free magnetic energy"). Nonlinearity might 
reflect m agnetic interactions a fter emergence (as would occur in an avalanche model 
of flares. ILu fc Hamiltonlll99ll ). or an aspect of the processes (dynamo, emgergence) 
that produce active region magnetic flelds. 



Given i?'s extensive scaling, it is interesting that ISchrijverl (120071 ) found that 
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ARs with similar values of R could produce flares with a range of energies. This 
suggests that the amount of flux near SPILs does not directly imply the amount of 
energy that will be released in any subsequent flares. It is possible that R scales 
with the total amount of magnetic energy available to be released, but that the size 
of the actual e nergy release d depe nds upon the details of the coronal magnetic field 



configuration. IWelsch fc Lil (120081 ) found that R tended to increase while total AR 
flux was increasing, implying that concentrations of flux near SPILs tend to arise 
by the emergence of flux, as opposed to the convergence of alr eady emerged flux. 



Emer gence has long been supposed to be a trigger of flares (e.g.. iPriest &: Heyvaerts 



19741 ). So, despite its extensive scaling, a large value of R might indicate the presence 
of a flare trigger, as well as the total amount of free magnetic energy present in a 
configuration. 

The flow property most strongly associated with average flare flux, Sr, has units 
of an energy flux: energy per unit area per unit time. In Table [5l raw values for this 
quantity are ~ 10^ km sec~^; mulitplying by both 10^, to convert km to cm, and 
the MDI disk-center pixel area, da ~ 2 x 10^^ cm^, and dividing by An (for consistency 
with the Poynting flux in equation [TS]) gives a power on the order of 10^^ erg sec~^. 
Integrated over the ~ 10^ seconds in a day or so, this represents the transport of 
~ 10^^ ergs day~^ of magnetic energy. While we have referred to Sr as a proxy 
for the integrated Poynting flux of magnetic energy across the photosphere, actual 
correlation between these quantities has not been demonstrated. If the true Poynting 
flux were an order of magnitude smaller, so ~ 10^^ ergs day~^, this would represent 
a enough ma gnetic energy to powe r a lar ge flare. This would contradict the view 



advocated by iMcClymont fc Fisherl (119891 ). who argued that even organized, near 



surface flows could not inject enough free magnetic energy into the coronal field to 
power large flares. Instead, they suggested that the magnetic energy that is released in 
very large flares must be carried into the corona with magnetic fields as they emerge. 

The view that flares are powered by free magnetic energy that emerges with 
AR fields is, however, supported by the fact that active regions which have recently 
flared are more likely flare again, compared to active regions tha t have not recently 



flared . This phenomenon has been referred to as persistence (jZirin &: Marquette 



199l[ ). The alternative view requires post-emergence energy buildup in the corona, 
and is plausible: flows could act on already emerged fields, leading to the storage of 
free magnetic energy in the corona, which could then be released in a flare. But why 
would flows preferentially inject magnetic energy into some ARs, making these flare 
more frequently, and not into other ARs? Indeed, it seems that some ARs emerge 
carrying a lot free magnetic energy, which is released in many flares; such ARs might 
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be "bom bad" fe.g.. lZirin fc Liggettlll987h . 



But if flare-prone ARs are born bad, why is average flare flux correlated with 
any photospheric flow property? It could be that Sr is in some way related to the 
emergence of new, non-potential magnetic fields. But if emergence were truly the 
dominant manner of free energy injection into the corona, then $ would be the more 
fundamental magnetic field property, and therefore presumably a better predictor of 
average flare flux. (Perhaps our timescale for computing $ — nominally, 96 minutes 
— is too short; see below for more on this point.) It could be that those ARs 
that emerge with the largest amount of free magnetic energy, due to the presence of 
electric currents, are also the most likely to exhibit large flows — perhaps as a result 
of Lorentz forces arising from these currents. 

In addition, it is possible that cancellation, driven by seemingly random flows, 
might lead to the buildup of free magnetic energy in the corona , which might then be 
releas ed in flares and CMEs. This process was not considered by lMcClymont fc Fisher 
(119891 ). Turbulent convective flows tend to systematically drive separate photospheric 
magnetic features int o close proximity in downflow lanes between supergranules (e.g., 
Schrijver et al.l 119971 ): while same-sign features simply merge, opposite-sign features 
undergo cancellation. Barring unrealistic symmetry, such cancellation generally forces 
reconnection between coronal fields, which has the effect of "stranding" any magnetic 
helicity that was present in those field s in the coro na, and tends to increase the 
magnetic free energy of the coronal field (I WelschI 120061 ). The buildup of free magnetic 
energy by this process, however, might be too slow to explain many large flares, which 
can recur frequently in the same active region. 

We were unable to clearly associate any intensive properties of flows with the 
occurrence of flares. Taken at face value, this might suggest t hat the flare mechanism 
itself is fundamentally extensive, as in the avalanche model of iLu fc Hamilton! (jl99ll ). 
for example. We must, however, qualify this result in several ways. As discussed 
above, our measure of flaring activity is extensive, which probably biases our anal- 
ysis against detecting relationships between intensive flow properties and flares. In 
addition, our flow estimates are almost certainly subject to serio us systematic errors . 
Tests of FLCT, DAVE, and other flow estimation techniques by lWelsch et al.l (120071 ) 
with sequences of synthetic magnetograms from MHD simulations, in which the true 
velocities are known, showed that both FLCT and DAVE systematically underesti- 
mated the fluxes of magnet ic helicity an d energy. (These shortcomings led to the 
development of DAVE4VM Schuckil2008| . which uses information about the horizon- 
tal magnetic field — unavailable in this study — to make improved estimates of flows 
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and the fluxes of magnetic energy and lielicity.) Moreover, tfiese tests also showed 
that FLCT systematically underestimated flow speeds. So the tracking methods and/ 
or data (2" spatial resolution, 96-min. cadence) employed in this study might be too 
limited to reveal any intensive flow patterns that are associated with flares. Or it 
might be that intensive flow properties are important, in some cases, for increasing 
coronal free magnetic energy, and/ or for triggering flares, but that the efficacy of 
a given flow pattern at injecting significant energy or triggering flares depends upon 
the details of the particular photospheric and coronal magnetic field configuration. 
In this case, statistical relationships between flows and flares would be hard to iden- 
tify. This last possibility motivates development of time-dependent models of coronal 
magnetic field evolution, driven by magnetograms, as tools necessary to understand 
(and eventually predict) the coronal response to photospheric flows. 

Could it be that some flow patterns are, in fact, anti- correlated with flare onset? 
Some previous results have suggested an inverse correlation between flow vorticity 
and flare occurrence. In his talk at the SHINE meeting in 2007, A. A. Pevtsov 
presented preliminary results suggesting that flares are more likely to occur after a 
time interval in which sunspots were observed to rotate. He presented a basic physical 
picture to explain these observations: (1) the solar interior acts as a current driver 
that imposes a fixed amount of twist (equivalently, current) on the flux systems that 
comprise active regions; (2) if the twist in the coronal portion of a flux system does 
not match the tw ist on the portion in the in terior, twist will propagate to equilibrate 



the two portions (ILongcope &: Welschll2000l ): (3) the twist propagation will be seen as 



a rotation of photospheric footpoints, which will slow as the twist gradient across the 
photosphere shrinks; (4) the twist injected into the corona implies that free magnetic 
energy is present in the corona, which will be released in a flare or flares. Based upon 
this model, one should expect that, statistically, flares will be more likely when the 
corona has more twist, which should occur when footpoint rotation has begun to slow 
or has ceased. 

This preliminary investigation of the statistical relationships between flows and 
flares has provided useful insights. We have investigated the properties of AR photo- 
spheric flows, characterized the correlation times of AR fields and flows, and identified 
a quantity, Sr, that might be useful for flare prediction. 

Further study will be necessary to fully understand the role of photospheric 
evolution in solar flares. A study similar to that conducted here, but with a larger 
sample of ARs, could verify that our results are not statistical aberrations, due to our 
small sample size. Also, our tracking techniques are primarily senstive to flows on the 
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spatial scale of the windowing parameter used by DAVE and FLCT, ~ 17". Flow-flare 
relationships on smaller and larger spatial scales should also be investigated; perhaps 
particular flows on very small scales, or coherent flows on large scales, are important. 

In addition, evolution on other time scales should also be investigated. Our 
correlation analysis suggested that flux emergence, quantified by the change in flux 
$ over our nominal 96-minute magnetogram cadence, was correlated with average 
flare flux, but not as strongly as other field and flow properties that we considered. 
But fractional changes of flux in magnetograms are not large on this time scale, 
which is short relative to the characteristic hours-long emergence timescale of active 
region flux. Consequently, the statistical relationship between flux emergence and 
flares on longer time scales (e.g., 6-48 hours) could reveal stronger correlations with 
flare activity. In addition, the decorrelation rate of magnetograms (see Figs. M - 
[To!) separated in time by intervals longer than our nominal cadence could provide a 
quantitative measure of AR magnetic evolution that might also be correlated with 
flares. Work on these efforts is underway. 
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A. Mercator Projections 

A map projection is a two-dimensional, planar representation of a two-dimensional, 
spherical surface. In solar physics, instruments typically record a two-dimensional im- 
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age of the Sun in some wavelength range as a plane-of-sky (POS) map of intensity 
variations. When the recorded emission originates in a relatively narrow atmospheric 
layer — as in a typical photospheric magnetogram, for instance — this procedure 
amounts to a projection from the imaged, spherical atmospheric layer on the ob- 
server's side of the Sun onto the image plane. From the viewpoint of scientific analy- 
sis, this projection poses many problems; for instance, foreshortening distorts spatial 
relationships between image features. 

To ameliorate this and other undesirable properties of the POS projection, the 
POS positions of image data points can be converted to latitude and longitude in some 
heliographic coordinate system. The origin of the Stonyhurst heliographic coordinate 
system (HCS), for instance, is at the i ntersection of the apparent solar equator and the 



central meridian as seen from Earth (IThompsonl 120061 ). This conversion from POS to 



HCS coordinates is a deprojection. Unfortunately, after transformation into longitude 
and latitude coordinates (0, 0) for a given hehographic coordinate system, data that 
were regularly spaced in the POS projection are typically not regularly sampled in 
(0, 6). (We use 6 for latitude, to distinguish it from the co-latitude 6 typically used in 
spherical polar coordinates.) Consequently, analyses of spatial relationships between 
image features are still problematic. 

Quantitative analyses of the spatial properties of irregularly sampled data lying 
on a spherical surface therefore require interpolation of these data onto a regular 
coordinate system. While interpolation onto spherical gridpoints uniformly spaced 
in (0, 9) is possible, both the distances and directions between nearby gridpoints in 
such systems vary in complex ways in different regions of the grid. Instead, we have 
opted to use the Mercator projection to convert the data from (0, 0) to a rectilinear 
coordinate system. Often this two-projection procedure — deprojecting from the 
POS coordinates in the image plane to the HCS, and then transforming the data again 
by interpolating onto spherical points corresponding to a uniform grid — is loosely 
referred to as "deprojecting" the data, even though it is probably best described 
as a reprojection. In the discussion that follows, we assume the HCS positions of 
the POS image points have already been computed, using, e.g., the IDL functions 
XY2LONLAT.PRO in SSWIDL's /ssw/gen/idl/solar/ . 

We chose the Mercator projection because it is conformal: locally, it does not 
distort shape. This means the Mercator projection should be well-suited for unbiased 
(in direction) tracking of features in image sequences, or any study of spatial variations 
in image structure, e.g., the magnitude and direction of intensity gradients. But, since 
any projection of a spherical surface onto a plane must distort either shape or scale 
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(if not both!), the pixel scale varies over the projected image. (Equivalently, this 
projection is non-authalic, or not area-preserving.) We have developed our Mercator 
projection procedure for tracking studies, on the assumption that correcting distortion 
of scale is easier than correcting directional bias in tracked displacements introduced 
by a non-conformal deprojection. 

The distortion of scale in global Mercator projections is a function of latitude 
alone. The horizontal Mercator coordinate, x, is mapped one-to-one with the helio- 
centric longitude, (p, and is independent of latitude. Because the physical distance 
between lines of constant longitude decreases as cos(^^) with increasing latitude, the 
projection's scale (the actual distance corresponding to a fixed distance in the repro- 
jected image) must decrease with increasing latitude. Displacements in the vertical 
Mercator coordinate, dy, corresponding to a fixed latitudinal displacement, d9, in- 
crease toward the poles. Quantitatively, y is related to latitude 9, via 

y = ln (^tan(^) + sec(^)) . (Al) 

As — > 7r/2, therefore, y diverges logarithmically. Differentially, dy — sec{9)d9, and 
the dy corresponding to an interval d9 also diverges toward the pole. To compensate 
for these distortions, small distances — e.g., tracked displacements, or the separations 
between pixel centers when computing gradients — must be multiplied by a factor of 
cosine latitude, cos(^), when converting Mercator distances to physical distances. 

It is desirable to preserve the spatial resolution in the original data set. Typically, 
this means the pixel width, As, in the projected data should match the pixel size at 
the sub-observation point (disk center) in the POS data. Hence, this fixes the spatial 
sampling rate of the original image that must be used to interpolate the POS data 
onto the Mercator grid. The number of necessary longitudinal samples for an interval 
of longitude A0 is the horizontal distance corresponding to this interval, = i?© A0, 
divided by the original pixel size. As. 

To preserve spatial resolution in the latitidunal direction is more complex. For 
a given interval of latitude, A9 = 9f — 9i, the interval Ay = y f — y^ is given by 

A, = lnf^^^^(^l±^^ . (A2) 
\^ tan(^i) + sec(^0 J 

The number of spatial samples necessary to preserve the initial pixel scale As is 
therefore the distance corresponding to the interval of latitude, Lg — RqAB, divided 
by As, and multiplied by the ratio Ay/A9. For A^ e {30°, 45°, 75°}, the ratios 
Ay/A9 are {1.05,1.12,2.02}, respectively. Use of too few samples in the vertical 
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direction would distort shapes in the image, and destroy the projection's conformal 
properties. 

The necessary interpolation is done by first computing the Delaunay triangula- 
tion of the irregularly sampled data in (latitude-longitude), then interpolating onto 
the latitude-longitude points corresponding to the regular Mercator grid of the output 
array. (This output grid is rectilinear in [latitude- longitude] , but not regularly spaced 
in latitude.) 

When deprojecting sequences of image data, we found that small changes in 
instrument pointing could lead to changes of -|-/- 1 pixel when computing array sizes 
from pointing information in the header files. Consequently, we fixed the output 
arrays for an entire image sequence, as opposed to determining these separately for 
each image. 

To verify that the algorithm worked properly, we projected latitudinal and lon- 
gitudinal bands onto a POS image of a sphere, and then reprojected the POS image 
backward. 

With the availabihty of full-disk vector magnetogram data from SOLIS/VSM, 
SDO/HMI, and other instruments, conformal projections could find wider application. 
To view data near the poles without distortion seen in typical global-scale Mercator 
projections (in which the line of tangency between the sphere and projection cylinder 
lies along the equator), local Mercator projections (in which the line of tangency runs 
through the center of the deprojected region) could be used, in a manner similar 
to the universal transverse Mercator (UTM) approach used in terrestrial maps, e.g., 
NATO topographical maps. 
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Table 1. Active Region Sample 



NOAA 

AR # 


Start Time 


Stop Time 


Avg. 

(/' 


GOES Flux^ Avg. GOES Flux^ , > 

m-2 day-i) (/, A^^ m-- day-i) 





7981 


1996-08-01, 04:19 


1996-08-04, 


22:24 


0.8 


0.0 


7986 


1996-08-28, 04:51 


1996-08-31, 


22:24 


0.4 


0.0 


7997 


1996-11-22, 13:03 


1996-11-23, 


17:39 


0.0 


0.0 


7999 


1996-11-24, 06:27 


1996-11-28, 


20:48 


8.7 


8.0 


8003 


1996-12-08, 01:36 


1996-12-11, 


12:48 


0.5 


0.0 


8004 


1996-12-15, 08:03 


1996-12-20, 


06:27 


0.1 


0.0 


8005 


1996-12-16, 14:24 


1996-12-21, 


12:51 


0.6 


0.5 


8006 


1996-12-17, 08:00 


1996-12-21, 


15:59 


0.0 


0.0 


8015 


1997-01-30, 03:15 


1997-02-03, 


22:24 


0.3 


0.0 


8020 


1997-03-10, 20:48 


1997-03-15, 


03:12 


0.0 


0.0 


8038 


1997-05-09, 12:52 


1997-05-14, 


00:04 


0.3 


0.3 


8040 


1997-05-18, 08:04 


1997-05-23, 


12:52 


3.1 


2.9 


8052 


1997-06-13, 14:24 


1997-06-17, 


22:24 


0.0 


0.0 


8076 


1997-08-28, 03:11 


1997-09-01, 


09:39 


2.5 


1.9 


8078 


1997-08-29, 01:39 


1997-09-01, 


22:24 


0.0 


0.0 


8081 


1997-09-05, 01:39 


1997-09-09, 


12:48 


0.0 


0.0 


8082 


1997-09-05, 01:39 


1997-09-09, 


23:59 


0.2 


0.0 


8083 


1997-09-06, 14:23 


1997-09-09, 


12:48 


3.7 


3.3 


8084 


1997-09-08, 09:39 


1997-09-13, 


12:48 


1.0 


0.6 


8085 


1997-09-12, 01:36 


1997-09-14, 


06:24 


1.4 


1.2 


8086 


1997-09-15, 14:27 


1997-09-20, 


12:51 


1.4 


1.2 


8088 


1997-09-24, 08:03 


1997-09-26, 


22:24 


45.4 


44.2 


8087 


1997-09-26, 17:36 


1997-09-29, 


03:12 


0.0 


0.0 


8090 


1997-10-05, 11:15 


1997-10-07, 


06:27 


0.0 


0.0 


8093 


1997-10-10, 17:35 


1997-10-15, 


12:51 


0.0 


0.0 


8096 


1997-10-16, 17:35 


1997-10-22, 


08:03 


0.0 


0.0 


8097 


1997-10-19, 19:11 


1997-10-24, 


16:00 


0.7 


0.7 


8100 


1997-10-31, 01:39 


1997-11-03, 


23:59 


88.6 


87.2 


8108 


1997-11-17, 22:23 


1997-11-22, 


06:23 


3.1 


1.6 


8123 


1997-12-17, 01:39 


1997-12-21, 


22:17 


0.1 


0.0 


8135 


1998-01-11, 17:36 


1998-01-15, 


20:48 


0.0 


0.0 


8156 


1998-02-13, 17:36 


1998-02-18, 


22:24 


3.9 


2.8 
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Table 1 — Continued 



NOAA 


Start Time 


Stop Time 


Avg. 


GOES Flux^ 


Avg. GOES Flux^ , > G^ 


AR# 














m~^ day^-*^) 


(/i W m"^ day"-*^) 


8158 


1998-02-15, 


01:36 


1998 


-02-18, 


22:24 




0.1 


0.0 


8179 


1998-03-13, 


03:12 


1998 


-03-17, 


22:24 




20.6 


19.7 


8210 


1998-04-28, 


17:39 


1998 


-05-03, 


16:00 




46.4 


45.2 


8214 


1998-05-02, 


17:36 


1998 


-05-05, 


16:00 




11.2 


10.7 


8218 


1998-05-10, 


03:15 


1998 


-05-14, 


22:24 




2.5 


1.5 


8222 


1998-05-15, 


17:36 


1998 


-05-19, 


03:15 




0.5 


0.5 


8220 


1998-05-15, 


14:24 


1998 


-05-19, 


09:39 




0.6 


0.6 


8227 


1998-05-31, 


01:39 


1998 


-06-02, 


12:51 




0.0 


0.0 


8230 


1998-06-02, 


11:15 


1998 


-06-06, 


22:24 




0.0 


0.0 


8243 


1998-06-17, 


01:36 


1998 


-06-21, 


12:48 




4.6 


4.1 


8375 


1998-11-02, 


03:11 


1998 


-11-06, 


22:24 




20.0 


19.4 


8377 


1998-11-07, 


03:12 


1998 


-11-09, 


19:15 




0.0 


0.0 


8392 


1998-11-26, 


01:36 


1998 


-11-29, 


12:48 




2.6 


2.6 


8393 


1998-11-28, 


01:36 


1998 


-12-01, 


22:24 




0.4 


0.4 



^ After lAbramenkol (120051 ) 

^"> G" denotes only flare with peak flux > GOES Gl.O level 
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Table 2. Magnetic Field and Flow Properties Calculated 



Property 


Scaling^ 


Related Physical Quantity 




Citation(s)^ 






E 


Quantify AR "magnetic size" 




Fisher et al. 1998: 
Leka & Barnes 2003 




E 


Flux emergence, cancellation 




Priest & Hevvaerts 1974: 
Linker et al. 2003 




I 


Flux emergence, cancellation 




?j 


1$!/$ 


I 


Unsigned flux changes 




n 


M[Br] 


I 


Field moments 




Leka & Barnes 2003 


M[\Bn\] 


I 


Unsigned field moments 






M[u] 


I 


Flow moments 






M[u\d^\/^u] 


I 


Flux-weighted flow moments 




— 




E 


Amount of motion in AR 




— 




E 


Unsigned fiux-weighted speeds 




— 


M[\r-V X u|] 


I 


Unsigned vorticity 




Brown et al. 2003 


A^ [ d<l>(r- V X u)|/<l> J 


I 


Flux-weighted unsigned vorticity 




" 


Er-(Vxu) 


E 


Net vorticity 




" 


Elf-(Vxu)| 


E 


Unsigned vorticity 




" 


^|rf<l>|(r- V X u) 


E 


Flux-weighted vorticity 




" 


J2\d^{r-V X u)| 


E 


Unsigned fiux-weighted vorticity 




" 


A^[V- u] 


I 


Converging/diverging flows 




Linker et al. 2003: 
Wane 2006 


A^[|V-u|] 


I 


Unsigned conv./div. flows 




" 


A^[|d$|(V-u)/<l>J 


I 


Flux-weighted conv./div. flows 




" 


A^[|d$(V-u)|/<l>J 


I 


Unsigned, flux-weighted conv./div. 


flows 




E(V ■ u) 


E 


conv./div. flows 






Eiv-u| 


E 


Unsigned conv./div. flows 






EM$l(v-u) 


E 


Flux- weighted, signed conv./div. flows 




EM$(v-u)| 


E 


Unsigned, flux-weighted conv./div. 


flows 




-DcOF 


E 


Convergence of AR flux 




Wane 2006 


^COF 


E 


Shearing of AR flux 








E 


Flux near PILs 




Falconer et al. 2003: 
Schriiver 2007 


R 


E 


Change in flux near PILs 




)) 


R/R 


I 


Change in flux near PILs 




5) 
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Table 2 — Continued 



Property 


Scaling*^ 


Related Physical Quantity 


Citation(s)'^ 




E 


Poynting flux proxy 








Flows along contours of Br 


Demoulin & Bereer 2003 






Flux-weighted contour flows 








Flows along gradients of Br 








Flux-weighted gradient flows 








Converging motions along PIL 


Linker et al. 2003 


E^UgBr 


E 




11 


M[WUcBr] 


I 


Shearing motions along PIL 


Antiochos et al. 1999 


E WUcBr 


E 


11 




Y.BI 


E 


Component of Photospheric Magnetic Energy 






E 


Kinetic energy 




E \Br\u' 


E 


Field-weighted kinetic energy 





^I = intensive, E = extensive 

'^References chosen are representative, not necessarily either exhaustive or initial. 



Table 3. Magnetic Flux Correlations 





R 




k 











0.69 ^ 


0.18 


b 


0.19 '^ 


0.91 '^ 


R 


0.53 ^ 




0.60 


b 


0.32 '^ 


0.84 '^ 


$ 


0.24 ^ 


0.54 ^ 






0.48 ^ 


0.40 ^ 


R 


0.21 ^ 


0.54 ^ 


0.48 


a 




0.17 '^ 




0.88 


0.61 ^ 


0.50 


a 


0.23 ^ 





^ linear (Pearson) correlation 
''rank-order (Spearman) correlation 



Table 4. Flare Frequencies and Skill Scores 



Wind. 


Code 


Varial)le(s) Considerod [Disc. Coeff.] 




PNF/F 


PF/NF 


PNF/NF 


SS 






JJAV £-1 


ri 


/ y 


Zo 1 




zo4o 


U. it) 








1 09 


91 4 


fi7 


zozo 


1Q 
u. ly 






P f1 1 Ql '^r, fn 771 /lRr,l\ fn ^91 


1 nfi 


91 
z±u 


70 


9'^99 


90 




r VjKj i 


c 

•J it 


1 nQ 


91 Q 
zio 


/D 


9Q1 R 
zOiD 


n 1 R 

U.iD 






<?r, fl 91 1 R n 1 fil 


1 07 


90Q 
zuy 


70 


9'^99 


90 
u.zu 






'in ri fi4l /? ri 111 '^lRr,l?/2 [-0 471 


1 07 


9nQ 
zuy 


uo 


9'?94 
zoz^ 


90 
u.zu 


DiN 


JJAV III 


f? 
r£ 


QO 

oy 


sn 


fi9 
OZ 


9c;o7 


-U.Uo 






R [9 nil fn 7fil 


40 


7Q 
/ y 


79 
( z 


9p;i 7 
ZO i ( 


08 

-U.Uo 






R \^ QSl '^T^ [n fifll /lRr,l\ [(1441 
-ft [i.yoj, [U.DUj, \\nji\i [u.^^j 


49 


77 


fi7 


9^199 
zozz 


O^i 

-U.UO 




r ij^^ ± 


Jri 


oy 


oU 


fi9 
uZ 


9c;97 
ZOZ ( 


n ns 

-U.Uo 






p [9 ni 1 Q Fn 71 1 
K [^.UiJ, Or [U. ( ij 


/1 1 


78 
/ o 


DD 




n n7 

-U.U ( 






R [1 QQl [n '^'?1 d) [n 471 


49 


77 


fi4 




04 




JJAV III 




9nn 
zyu 


OOO 


Do 


1 0(^7 

iyb/ 


n QQ 
U.oo 






Q Fl 1 ol /I R l\ ri nal 

OR [i-iyj, \|-Di?|/ [i.usj 


o4o 


OZO 


1 ni 

iUi 


1 OQ/l 

iyo4 


n Q7 






5r [1.63], {\Br\) [0.95], Y,\Br\u^ [-0.39] 


349 


324 


90 


1945 


0.37 


24C 


FLCT 


Sr 


299 


374 


72 


1963 


0.33 






Sr [1.22], {\Br\) [1.03] 


335 


338 


76 


1959 


0.37 






[1.19], {\Br\) [0.71], ^[5^1] [0.37] 


330 


343 


76 


1959 


0.37 


24N 


DAVE 




107 


221 


83 


2297 


0.15 






Sr [1.19], i? [0.85] 


106 


222 


72 


2308 


0.17 






[0.82], R [0.78], [0.55] 


105 


223 


73 


2307 


0.18 


24N 


FLCT 


Sr 


107 


221 


74 


2306 


0.16 






Sr [1.34], [0.74] 


112 


216 


67 


2313 


0.18 






>Sr [1.35], i? [0.47], J^24,current [0.45] 


115 


213 


63 


2317 


0.20 



^PF = predicted to flare; F = flare(s) occurred; PNF = predicted no flare; NF = no flare occurred 



Table 5. Means of Selected Raw and Standardized Variables for Flaring and Non-Flaring Populations 



Window 


Variable [Units] 


, raw 


NF^ , raw 


F + NF, raw 


Std. Dev.-i 


F, standardized^ 


NF, standardized 


6C 


Sr [km G7s] 


3.6 xlO^ 


1.4 xlO^ 


1.6 xlO^ 


1.6 xlO^ 


1.3 


-1.7 xlO-^ 




R[G] 


7.3 xlO=^ 


1.0 xlO^ 


1.8 xlO^ 


4.5 xlO^ 


1.2 


-1.6 xlO-^ 


6N 


R[G] 


1.0 xlO^ 


1.4 xlO^ 


1.8 xlO^ 


4.5 xlO^ 


1.9 


-8.8 xlO-2 




Sr [km GVs] 


4.0 xlO^ 


1.5 xlO^ 


1.6 xlO^ 


1.6 xlO^ 


1.5 


-7.0 xlO-2 


24C 


Sr [km GVs] 


3.2 xlO^ 


1.1 xlO^ 


1.6 xlO^ 


1.6 xlO^ 


0.99 


-3.3 xlO-^ 




( \Br\ ) [G] 


5.1 xlQi 


3.1 xlQi 


3.6 xlO^ 


1.6 xlO^ 


0.97 


-3.2 xlO-^ 


24N 


Sr [km GVs] 


3.6 xlO^ 


1.4 xlO^ 


1.6 xlO^ 


1.6 xlO^ 


1.2 


-1.7 xlO"^ 




R[G] 


6.5 xlO^ 


1.1 xlO^ 


1.8 xlO^ 


4.5 xlO^ 


1.1 


-1.5 xlO"^ 



*FLCT flow estimates were used. 

'^Overlines denote averages; F = flare(s) occurred 

'^NF = no flare occurred 

■^Standard deviation of combined F-I-NF populations. 

''Derived by subtracting the mean from raw data, then dividing by the standard deviation. 
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Fig. 1. — FLCT flows overlaid on a grayscale image of Bji from AR 8038; white and 
black vectors represent flows in negative and positive regions, respectively. White and 

black contours of are ovcrplottcd in negative and positive regions, respectively, at 
100 G intervals between 100 G and 500 G. The time shown at top corresponds to the 
end of the 96-min. interval over which flows were estimated. 
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Fig. 2. — DAVE flows overlaid on a grayscale image of Bn from AR 8038; white and 
black vectors represent flows in negative and positive regions, respectively. White and 
black contours of Br are overplotted in negative and positive regions, respectively, at 
100 G intervals between 100 G and 500 G. The time shown at top corresponds to the 
end of the 96-min. interval over which flows were estimated. Note that the vector 
length scale differs from that in Figure [H 



Ijj 



RGn{(-Or:cje^ Corr: 0.55 



■■■■■■ 



-0,10-0.05 0.00 0.05 0,10. 



□near Ccrr: Cs. 47 
: Sor^k-Or-der Core: 0.56 

' I ■ ■ ■ 



O.-O 0.3 
FLtjT Li^ [km sec' 'J 



Fig. 3. — FLCT and DAVE flows from Figures [T] and [21 respectively, are compared in 
scatter plots of (left) and Uy (right). We also show linear and rank-order correlation 
coefficients between the flow components. 
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Fig. 4. — Examples of flow properties from FLCT flows corresponding to the flow field 
shown in Figure [TJ Top row: Grayscale maps of (left) uBj^ (right). Black and 

white contours of Br are overplotted in negative and positive regions, respectively, at 
100 G intervals between 100 G and 500 G, for all panels. Bottom row: Grayscale maps 
of I (i$ I -weighted versions of the horizontal divergence (left) and unsigned vertical curl 
(right). 




Fig. 5. — Top: The background grayscale shows Br (white is positive, black negative), 
and white and black contours of Br are overplotted in negative and positive regions, 
respectively, at 100 G intervals between 100 G and 500 G. Components of FLCT 
flows from Fig. [1] that lie along gradients and contours of Br are shown with blue 
and red vectors, respectively. For clarity, not every vector is plotted. Bottom left: A 
grayscale map of Wpn- and 5/j-weighted gradient flows, used to quantify converging 
motions. Black and white contours of Br are overplotted in negative and positive 
regions, respectively, at 100 G intervals between 100 G and 500 G. Red contours of 
the PIL- weighting map Wpn are overlaid. Bottom right: A grayscale map of Wpn- 
and i?/j-weighted contour flows, used to quantify shearing motions. All other image 
components are as with the bottom left image. 



■SKt-' new i=S!Sf) ■i.sem 
i^j fcrT' 




Fig. 6. — Combined log-scaled shading histogram and scatter plot of speed estimates 
from FLCT (top) and DAVE (bottom) as a function of pixel-averaged field strength 
Bji (see text). The grayscale for FLCT ranges from 2.4 - 6.2, and for DAVE from 2.4 
- 6.0. Higher speeds tend to be seen in weaker fields. 
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Fig. 7. — FLCT and DAVE rank-order correlation coefficients between whole-disk- 
passage averaged flow properties and flare power, averaged over the 46 active regions 
in our study. To the extent that correlation coefficients lie on the slope = 1 dashed 
line {not a fit), the averaged fiow properties from the FLCT and DAVE methods 
agree. The quantities correlated are numbered by distance from (0,0), farthest first; 
hence, smaller numerical labels imply stronger correlation. Many of the strongest 
correlations do not require estimating fiows. Even ranks are plotted above their 
symbols, and odd ranks below. 
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Fig. 8. — FLCT and DAVE linear correlation coefficients between whole-disk-passage 
averaged flow properties and flare power, averaged over the 46 active regions in our 
study. To the extent that correlation coefficients lie on the slope = 1 dashed line {not 
a fit), the averaged flow properties from FLCT and DAVE methods agree. The corre- 
lations are numbered by distance from (0,0), farthest ffist; hence, smaller numerical 
labels imply stronger correlation. As in rank-order correlations in Figure [H many 
of the strongest correlations do not require estimating flows. Even ranks are plotted 
above their symbols, and odd ranks below. 
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Fig. 9. — For ARs 7981, 8038, and 8210, frame-to- frame correlations in Bn with time 
are shown in black; the thick lines represent linear correlation coefficients between the 
current and previous frames, while the thin lines show linear correlation coefficients 
between the current and initial frames. Analogous correlations are also shown for 
FLCT's estimated & Uy in blue & red (resp.), with thick & thin lines denoting 
current-previous & current-initial correlations (resp.). In addition, flare start times 
are shown by letters representing GOES flare class; the vertical placements of letters 
do not scale with actual flare fluxes. 
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Fig. 10. — As with Figure [9], but with correlations using DAVE's & Uy in blue & 
red, resp. 
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Fig. 11. — Top row: Scatter plots of rank-order (R.O.) and linear correlation coef- 
ficients (left and right panels, resp.) between average flare power and magnetogram 
or flow properties from FLCT (horizontal axis) and DAVE (vertical axis), within the 
6-hour time window, 6C, centered on the average time of the magnetograms used to 
estimate flows. Quantities are labeled in order of decreasing distance from (0,0), the 
point on the plot which corresponds to complete lack of correlation with average flare 
power. Low-ranked properties are not shown. Odd rank numbers appear below their 
corresponding plot symbols, while even rank numbers appear above their symbols. 
To the extent that flare correlations with flow properties from FLCT and DAVE flow 
agree, they lie along the dotted hne (not a fit). Quantities that do not require a flow 
estimate (e.g., total fliix, $) lie on the diagonal line. Bottom row: same as for top, 
but for the 24-hour time window 24C. Some of the quantities most strongly correlated 
with average flare power — e.g., R — do not require estimating flows. Our proxy for 
the Poynting flux, 5*^, exhibits comparable correlations. 
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Fig. 12. — As with Figure fTTl these plots show correlation coefficients between proper- 
ties of magnetograms or flows from FLCT (horizontal axis) and DAVE (vertical axis), 
and average flare power from the 6N (top row) and 24N (bottom row) windows. In 
each panel, quantities are labeled in order of their distance from (0,0). Low-ranked 
properties are not shown. Odd rank numbers appear below their corresponding plot 
symbols, while even rank numbers appear above their symbols. As with Figure [TTl 
many quantities strongly correlated with average flare flux do not require estimating 
flows, but our proxy for the Poynting flux, Sr, exhibits comparable correlations. 
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Fig. 13. — Squares show log^Q average flux from flares above GOES Cl.O level while 
each AR was within 45° of disk center vs. log average total unsigned flux in these 
ARs over that time interval. Horizontal bars show log of RMS variations flux during 
that interval; errors from magnetograph noise for any single flux measurement are 
about the size of the plotted symbols or smaller. Triangles show log of magnetic flux 
in the 22 ARs that did not flare above GOES Cl.O level. We also list: the hnear and 
rank-order (R.O.) correlation coefficients between the raw magnetic and flare fluxes 
for all ARs; the hnear and R.O. coefficients between the logs of magnetic and flare 
fluxes for flaring ARs; and the slope of a least- absolute- deviation fit (dotted line) 
of the log of flare flux as a function of log of magnetic flux for ARs that flared above 
the Cl.O level. 
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